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ABSTRACT 


In  an  effort  to  reduce  the  time  and  cost  of  ballistic  experimental 
development  projects,  the  factorial  experiment  and  Box  technique  have 
been  applied  to  ballistic  experiments  with  aircraft  seat  ejection  cata¬ 
pults  and  rockets,  high-low  guns,  Davis  guns,  and  recoilless  rifles. 

This  has  resulted  in  a  reduction  in  the  number  of  rounds  fired,  with 
little  or  no  reduction  in  the  validity  of  the  analysis  of  variance. 

A  detailed  discussion  of  application  of  the  Box  technique  to  the 
factorial  data  is  presented.  This  application  results  in  the  determina¬ 
tion  of  a  "zone  of  suitable  performance"  which  makes  use  of  interaction 
effects  to  provide  greater  flexibility  in  the  selection  of  design  parameters. 
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APPLICATION  OF  FACTORIAL  EXPERIMENT  AND  "BOX"  TECHNIQUE 

TO  BALLISTIC  DEVICES 


INTRODUCTION 


In  the  experimental  development  of  ballistic  systems  at  Frankford 
Arsenal  we  are  faced  with  a  wide  variety  of  experimental  problems. 

For  example,  in  recent  years  we  have  been  concerned  with  recoilless 
weapon  systems,  aircraft  seat  ejection  catapults  and  rockets,  thrusters, 
high- low  guns,  and  reactionless  launchers.  Some  of  these  systems  are 
required  to  function  repeatedly  with  performance  variations  of  the  order 
0.  1  percent;  others  are  one-shot  devices  which  must  function  reliably 
with  performance  over,  under,  or  within  certain  prescribed  limits. 

Sample  size  varies  a  great  deal,  as  well  as  the  type  of  performance 
requirement.  In  development  of  items  such  as  small  arms  cartridges 
we  can  fire  thousands  of  experimental  rounds,  while  a  seat  ejection 
catapult,  for  example,  limits  us  to  20  or  30  to  50  rounds. 

By  use  of  factorial  experimental  design  techniques  and  analysis, 
combined  with  physical  interpretation  of  the  data  in  terms  of  response 
surfaces  {the  "Box"  technique)  as  suggested  by  Dr.  Box,  ^  a  treme:\dous 
flexibility  of  standard  statistical  practices  has  been  achieved.  This 
method  has  been  applied  in  one  way  or  another  to  the  devices  mentioj.’ed 
previously.  As  examples,  our  studies  with  the  reactionless  launcher, 
an  analog  computer  simulation  of  a  thruster,  and  the  "box"  of  a  seat 
ejection  catapult  will  be  discussed. 

The  presentation  herein  illustrates,  in  chronological  order,  a  step 
by  step  experimental  evaluation  of  the  technique.  (The  experimental 
evaluation  was  preceded  by  an  abstract  evaluation  which  is  not  reported 
here.)  First,  existing  data  from  a  seat  ejection  catapult  development 
was  studied  to  determine,  in  a  prelimmary  way,  the  method's  effective¬ 
ness,  the  required  type  of  experiment,  and  some  of  the  experimental 
pitfalls.  Second,  a  report  is  made  of  the  theoretical  study  of  a  thruster, 
from  which  we  learned  something  about  the  response  surfaces  and  methods 
of  interpolation.  Finally,  a  discussion  is  presented  of  the  reactionless 
launcher  study.  This  latter  study  was  conducted  from  start  to  finish  using 
the  experimenta.l  design  methods  we  propose. 


*Box,  G . E .  P.  ,~^The  Exploration  and  Explanation  of  Response  Surfaces: 
Some  General  Considerations  a.nd  Examples Biometrics,  Vol  10, 

No.  1,  Mar  1954 
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APPLICATION  OF  THE  "BOX"  TECHNIQUE  TO 
M5(MOD)  CATAPULT  DATA 


The  possibility  of  applying  the  Box  technique  to  existing  data  for 
the  modified  M5  seat  ejection  catapult  was  considered.  Although  a 
carefully  controlled  experiment,  as  performed  in  the  reactionless 
launcher  study  (to  be  discussed  later),  is  required  to  obtain  fully  valid 
results,  a  preliminary  analysis  of  existing  data  by  the  Box  technique 
was  expected  to  give  some  indication  of  its  effectiveness.  Data  from 
24  firings  of  the  modified  M5  catapult  were  analyzed  using  three  vari¬ 
ables:  temperature  (T),  charge  (C),  and  web  (W),  each  at  two  levels 
for  two  propellant  compositions  (Lots  5655,  1  and  5656,1), 

At  that  time  the  modified  M5  catapult  was  to  meet  the  following 
requirements: 

1,  The  peak  acceleration  (g)  and  the  rate  of  change  of  acceleration 
(^)  were  not  to  exceed  25  g's  and  300  g/second,  respectively; 

2,  The  final  velocity  (v)  was  to  equal  or  exceed  80  fps. 

The  least  square  method  was  employed  to  fit  plane  surfaces*!*  to 
the  experimental  data  for  g,  g,  and  v,  yielding  the  following  equations: 

g  =  -308,  3W  +  0. 13C  +  0.097T  +  46.  1 

g  =  -1354W  +  0.358C  +  1,463T  +  287.6 

v=  152,  IW  +  0.3075C  +  0.0908T  +  60.08 

where  W  is  in  inches;  T,  in  °F;  and  C,  in  grams. 

The  equations  were  plotted  for  constant  values  of  C,  W,  and  T; 
i,e,,  the  intersections  of  the  g,  g,  and  v  responses  with  the  six  planes 
formed  by  choosing  constant  values  of  C,  W,  aildT  were  graphed 
(Figures  1  through  3).  The  lines  on  these  graphs  represent  the  inter¬ 
section  of  the  re:^ponse  surface  with  the  constant  planes.  For  example, 
Figure  3A  depicts  the  intersection  of  the  g,  g,  and  v  response  surfaces 
with  the  plane  formed  by  taking  the  temperature  as  70 °F.  The  arrows 
indicate  the  direction  of  increasing  magnitude  of  the  v  response  surface 
and  decreasing  magnitude  of  the  g  and  g  surfaces. 

*!'The  functions  are  not  really  plane  surfaces.  To  simplify  the  calculations, 
a  limited  range  of  the  parameter  is  chosen  so  that  the  variables  can  be 
considered  a  linear  function  of  the  parameters  within  that  range.  Caution 
must  therefore  be  exercised  when  interpolating  or  extrapolating.  For 
example,  the  origin  (W=T=:C~0)  is  not  a  valid  point  on  these  planes. 
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Constant  Response  Surface 
A  -  Charge  =  1?0  gm 
B  -  Charge  =  130  gm 
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figure  1 


Charge  1lel|;ht  Cs^« ) 
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Figure  3.  Constant  Response  Surface 
A  -  Temperature  =  70“  F 
B  -  Temperature  =  160“  F 
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The  next  step  was  to  form  the  six  constant  planes  into  a  box.  The 
response  surfaces  within  the  cube  were  obtained  by  joining  the  corre¬ 
sponding  curves  for  g,  g,  and  v.  Photo  1  (Appendix)  shows  this  box. 

The  thickness  of  the  response  surfaces  is  a  result  of  round  to  round 
variation  in  ballistic  performance.  This  illustration  is  qualitative; 
actual  thickness  must  be  determined  from  analysis  of  variance  of  the 
data. 

An  operating  point  (W-  ,  C^,  T^ )  which  satisfies  performance 
requirements  for  this  model  must  be  within  the  cube  volume  defined  by 
the  three  response  surfaces.  It  is  seen  that  the  g  and  g  requirements* 
are  not  met  by  all  points  within  this  space,  with  the  exception  of  points 
in  front  of  these  planes  (in  the  direction  of  the  arrows).  For  example, 
the  coordinates  of  point  Wq  =  0. 158  in, ,  Cq  =121  gm,  and  Tq  =  85®  F , 
give  a  web,  charge,  and  temperature  at  which  acceleration  change  is 
less  than  300  g/sec,  with  a  velocity  greater  than  80  fps.  We  see  further 
that  there  is  a  volume  surrounding  this  point  over  which  the  specifica¬ 
tions  will  be  met.  We  will  call  this  volume  the  zone  of  suitable  response. 
It  has  limiting  values  determined  by  the  geometry  of  the  response  sur¬ 
faces, 

A  better  operating  point  might  be  found  by  extending  the  v,  g,  and 
g  response  surfaces  outside  the  limits  of  the  box.  For  example,  it 
appears  that  a  new  constant  web  plane  for  webs  greater  than  W  =  0.16 
inch  will  increase  the  temperature  range  over  which  desired  perform¬ 
ance  is  achieved. 

In  addition,  the  response  surfaces  may  be  extended  in  the  direction 
of  increasing  or  decreasing  charge  or  temperature;  thus,  a  volume 
space  can  be  obtained  over  any  desired  range  of  web,  charge,  and 
temperature  on  the  basis  of  a  relatively  few  firings.  (Other  values, 
such  as  internal  volume,  expansion  ratio,  etc.,  could  be  used  instead 
of  those  chosen  for  this  particular  model.)  Any  extension  of  the  re¬ 
sponse  surfaces  outside  the  cube  which  represents  experimental  values 
is  only  as  valid  as  the  assumption  that  the  response  surfaces  are  planes. 
It  becomes  important,  then,  to  learn  something  about  the  response  sur¬ 
face.  In  particular,  the  hazards  involved  in  interpolation  and  extrapola¬ 
tion  should  be  studied,  A  start  was  made  in  this  direction  with  a  theo¬ 
retical  study  of  a  thruster. 


- - - - 2 - 

='<g  <  25  ft/second  and  g 


S300  g/second 
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THEORETICAL  STUDY  OF  A  THRUSTER 


An  analog  computer  was  used  to  develop  theoretical  response 
surfaces  for  a  thruster  which  moves  a  500-lb  load  vertically » Two 
restrictions  were  Imposed; 

1.  Maximum  pressure  was  to  be  less  than  7000  psi; 

2.  Final  velocity  was  to  be  greater  than  7,5  fps . 

About  60  computer  runs  were  made  for  various  design  parameter 
Combinations.  The  ballistic  design  parameters  which  were  considered 
were:  charge  (C),  propellant  web  (W)»  and  chamber  volume  (Vj,  ),  The 
intersection  lines  of  the  response  surfaces  with  the  planes  were  obtained 
graphically  from  the  results  of  the  60  simulations. 

Figure  4A  illustrates  the  intersection  of  the  response  surfaces  with 
the  plane  C  =  3  grams,  while  Figure  4B  is  the  intersection  with  the  plane 
=1,3  in.^  and  Figure  4C,,  the  intersection  with  the  .plane  W  ~  0,11  in. 

The  three  dimensional  representation  of  the  two  response  surfaces 
(pressure  =  7000  psi  and  velocity  =7,5  fps)  are  shown  in  Photo  2  (Ap¬ 
pendix),  Some  warping  of  the  response  surfaces  can  be  seen.  This 
illustrates  a  nonlinear  response.  However,  the  nonlinearity  is  well 
behaved.  No  oscillations,  peaks,  or  humps  occur.  A  linear  interpola¬ 
tion  should,  therefore,  be  adequate  if  the  box  is  small  enough.  At  most, 
second  order  terms  would  be  necessary.  The  size  of  the  box  should  be 
small  compared  to  nonlinearities,  but  large  compared  to  nonuniformities . 

Preliminary  experi.mental  work  in  ballistic  development  should  be 
directed  toward  determining  linearity  and  uniformity.  This  information 
is  essential  before  iCtting  up  the  factorial  experiment  so  that  the  dif¬ 
ferences  in  performance  levels  will  be  significant  and  so  that  the  com¬ 
plexities  of  nonlinear  interpolation  of  the  data  can  be  avoided.  In  addi  ¬ 
tion,  this  information  should  give  some  idea  of  the  range  of  validity  of 
extrapolations.  However,  it  is  a  good  practice  always  to  verify  extra¬ 
polation  experimentally.  (Proper  preliminary  work  should  eliminate  the 
need  for  extrapolation, ) 


“^'Details  of  computer  simulation  of  ballistic  devices  can  be  found  in  the 
following  references: 

Frankford  Arsenal  Report  R-1313,  "An  Analog  Computer  Study  of 
Interior  Ballistics  Equations,"  March  1956,  by  I;.  Stout  and  W.  A. 
Dittrich 

Frankford  Arsenal  Report  M61-17-1,  "Analog  Computer  Study  of 
Interior  Ballistics  of  Propellant  Actuated  Devices,"  April  1960,  by 
R.  Boritz  and  L.  Narisi 
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The  operating  volume  of  zone  of  suitable  response  is  seen  to  be 
triangular  in  cross  section,  opening  up  in  the  direction  of  increased 
chamber  volume  and  corresponding  increased  web.  Thus,  for  an 
increased  chamber  volume,  the  range  of  web  and  charge  over  which 
the  two  restrictions  would  be  met  is  greater.  Picking  a  set  of  values 
for  C,  W,  and  Vc  approximately  in  the  center  of  the  zone  of  suitable 
performance  would  thus  minimize  the  chance  of  violating  our  restric¬ 
tions  because  of  performance  variations  that  result  from  manufacturing 
tolerances.  A  larger  chamber  volume  would  allow  substantial  reduc¬ 
tion  of  these  tolerances.  The  actual  chamber  volume  allowable,  of 
course,  is  subject  to  the  physical  size  of  the  thruster  and  other  ballistic 
considerations,  such  as  ignition  and  expansion  ratio. 
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REACTIONLESS  LAUNCHER  STUDY 


The  reactionless  launcher  is  a  Davis  type  recoilless  gun  for 
ejecting  masses  from  a  ballistic  missile  during  flight.  The  launcher 
holds  two  projectiles,  as  shown  in  Figure  5.  In  the  particular  pro¬ 
ject  to  be  discussed  here,  the  masses  were  intended  to  decoy  anti¬ 
missile  missiles. 

The  decoys  are  of  many  sizes  and  weights,  and  are  launched  at 
a  wide  range  of  velocities.  The  weight  range  considered  was  20  to 
60  pounds,  and  the  velocity  ranged  from  50  to  110  fps.  The  wide 
range  of  performance  required  two  types  of  interior  ballistic  systems, 
direct  and  high-low  (shown  in  Figure  6),  There  are  two  types  of  pro¬ 
jectile;  the  bullet  type  (full  caliber),  which  fits  directly  in  the  bore  of 
the  gun,  and  the  spigot  type,  which  has  a  rod  that  fits  in  the  gun  barrel, 
with  the  payload  outside  the  gun.  These  are  shown  in  Figure  7. 

The  entire  study  involved  a  total  of  eight  variables:  ballistic  sys¬ 
tem  geometry,  decoy  geometry,  charge  weight,  decoy  weight,  shot- 
start  static  breaking  pressure,  expansion  ratio,  propellant  web,  and 
orifice  area  (high-low  system). 

To  blindly  set  up  a  factorial  experiment  at  two  levels  would  re¬ 
quire  the  firing  of  2°,  or  256,  rounds.  Replicating  three  times,  which 
is  reasonable  for  this  type  of  study,  would  lead  to  firing  more  than  750 
(of  the  order  10^)  rounds.  Instead,  we  isolated  factors  with  no  inter¬ 
actions,  such  as  the  type  of  chamber  -  the  high-low  chamber  was 
studied  separately  from  the  direct  chamber.  We  divorced  the  spigot 
projectile  from  the  bullet  type  for  the  direct  system,  but  not  for  the 
high-low  system  since  the  high-low  performance  would  not  be  ex¬ 
pected  to  depend  strongly  on  the  type  of  projectile.  As  a  result  of  this 
disjoining  process,  the  study  was  split  into  three  programs:  A,  B, 
and  C« 

In  program  A,  a  high-low  chamber  was  used  with  a  bullet  type 
projectile.  The  main  variables  were:  charge  weight,  shot- start 
static  breaking  pressure,*  and  orifice  area  of  the  high  pressure  cham¬ 
ber . 


In  program  B,  a  direct  chamber  was  used  with  a  bullet  type  pro¬ 
jectile  .  The  main  variables  were:  charge  weight,  propellant  web, 
and  shot-start  static  breaking  pressure. 


*Shot-start  is  a  rod  which  restrains  projectile  motion  until  chamber 
pressure  reached  a  predetermined  level. 


10 


36.231.S0367/ORD.60 


Ml  Vi  =  Mz  V2 

NET  REAaiON  » 


0  I 


lU 


36.231.S0364/ORD.60 


Figure  7,  Launchers,  Full  Caliber  and  Spigot  Type 


In  program  C,  a  direct  chamber  and  spigot  projectile  were  used. 
The  main  variables  were:  decoy  weight,  shot-start  static  breaking 
pressure,  and  spigot  design  (i.e.,  expansion  ratio). 

We  fired  factorial  experiments  at  two  levels  for  these  variables 
(eight  rounds  for  each  program).  For  the  three  programs  (A,  B,  and 
C),  which  were  replicated  three  times,  we  fired  a  total  of  8  x  3  x  3, 
or  72  rounds,  a  reduction  by  a  factor  of  10  in  the  number  of  rounds 
required.  The  change  in  the  magnitude  of  t  in  the  "Student  t  Test"  is 
less  than  5  percent  for  this  change  in  the  number  of  degrees  of  freedom 
and  for  confidence  levels  as  high  as  99.95  percent. 

The  discussion  is  being  confined  to  the  C  program,  as  this  amply 
illustrates  the  important  points  and  the  other  programs  are  similar. 

The  statistical  method  used  is  found  in  Kempthorne,*  The  data 
taken  were  peak  chamber  pressures,  peak  acceleration,  and  the  muz¬ 
zle  velocities  of  the  projectiles.  In  addition,  several  other  ballistic 
parameters  (such  as  piezometric  efficiency  and  ballistic  efficiency) 
were  examined.  Each  result  was  treated  separately,  in  the  manner 


*Kempthorne,  Design~and  Analysis  of  Experiments,  New  York;  John 
Wiley  &(  Sons,  IncT,  1^52 
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outlined  in  Kempthorne,  to  obtain  the  effects  of  each  variable  and  the 
interactions  between  variables.  The  significance  of  these  values  was 
ascertained  by  the  use  of  the  standard  error  and  "t”  test  at  both  the 
5  and  1  percent  levels. 

The  values  of  the  variables  investigated  in  this  program  are  shown 
in  Table  I;  test  results  obtained  are  shown  in  Table  II.  Results  of  the 
factorial  analysis  are  presented  in  Tables  III,  IV,  and  V,  showing  the 
effects  and  interactions  of  the  variables  on  peak  pressure,  peak  acceler¬ 
ation,  and  muzzle  velocity,  respectively. 

In  Tables  III,  IV,  and  V,  a  capital  letter  is  used  to  represent  the 
average  effect  of  the  corresponding  parameters.  In  Table  III,  for 
example,  P  =  -2490  psi  represents  the  difference  between  the  average 
peak  pressure  of  aH  rounds  fired  with  a  closed  spigot  (closed  spigot 
indicates  large  expansion  ratio;  consequently,  this  was  considered  the 
upper  level  of  this  parameter)  and  all  rounds  fired  with  an  open  spigot. 
Two  capital  letters  written  together'T'Wj.  P,  for  example)  represent  the 
interactions  of  the  two  cor  re 
Table  III,  P  =  -995  psi;  i. 
and  interactions  follows. 

The  interpretation  of  effects  and  interactions  is  as  follows.  The 
main  effect  (P,  for  example)  is  the  effect  on  the  variable  (pressure,  in 
Table  III;  acceleration,  in  Table  IV;  and  velocity,  in  Table  V)  of  in¬ 
creasing  expansion  ratio  (changing  from  closed  spigot  ^  open  spigot) 
averaged  over  all  possible  combinations  of  projectile  weight  and  shot- 
start  values.  It  is  desirable  now  to  determine  the  effect  of  expansion 
ratio  averaged  over  all  shot-start  values,  but  at  the  low  projectile 
weight.  This  is  denoted  symbolically  by  P  -  PW^  . 

In  Table  III,  for  example,  P  -  PW^  =  -2180  psi  indicates  that  using 
data  for  20-lb  projectile  weight  only  and  averaging  over  aU  shot-start 
values,  the  peak  pressure  is  reduced  2180  psi  in  changing  from  large 
expansion  ratio  (closed  spigot)  to  small  expansion  ratio  (open  spigot). 

For  data  from  the  60-lb  projectile  weight  and  all  shot-start  values 
(symbolically,  P  +  PWj),  we  have  -2790  psi.  The  fact  that  P  -  PW^ 
differs  from  P  +  PW^  indicates  an  interaction  between  projectile 
weight  and  expansion  ratio. 

The  results  in  Table  III  show  that  P  -■  200  psi  and  - 

W(.  P  =  810  psi.  Therefore,  the  projectile  weight  effect  when  the  open 
spigot  is  used  is  200  psi;  when  used  with  the  closed  spigot,  the  pro¬ 
jectile  weight  effect  is  810  psi.  The  difference  value  of  610  psi  (2790 
psi  -  2180  psi  and  810  psi  -  200  psi)  is  the  interaction  effect  between 
expansion  ratio  and  projectile  weight. 


spending  parameters.  Using  data  from 
f..,  ~  ^490  ,  Interpretation  of  effects 
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For  a  pictorial  representation  of  the  results,  the  variables  are 
laid  out  as  the  axis  of  a  transparent  cube.  The  corners  of  the  cube 
represent  the  eight  combinations  of  variables  fired.  The  yields 
(velocity,  acceleration,  and  peak  chamber  pressure)  are  assumed 
to  vary  along  the  edges  of  the  cube  according  to  the  predictions  of 
ballistic  theory.  Thus,  the  yields  at  the  corners  are  interpolated 
to  obtain  planes  of  constant  response.  (Ideally,  an  analog  computer 
analysis  is  desirable  to  calculate  the  planes  exactly,  as  was  done  for 
the  thruster  previously  discussed.) 

The  planes  indicated  in  Photo  3  (Appendix)  represent:  peak  pres¬ 
sure,  2800  psi;  velocity,  108  fps;  and  peak  acceleration,  360  g's. 
Points  within  the  transparent  cube  above  the  red  surface  (designated 
P)  represent  variables  which  result  in  pressures  below  2800  psi, 
Similarly,  points  in  front  of  the  V  surface  (green)  are  below  108  fps, 
and  behind  the  G  surface  are  less  than  360  g's.  Thus,  the  three  sur¬ 
faces  enclose  a  polygon  of  triangular  cross  section  which  is  the  zone 
of  suitable  response. 

Combinations  of  variables  near  the  surface  of  the  zone  may  re¬ 
sult  in  unsuitable  performance  as  a  result  of  round-to-rovind  varia¬ 
tions.  Analysis  of  variance  from  the  results  of  the  factorial  analysis 
and  interpolation  of  the  variance  along  the  cube  edge,  using  the  same 
technique  as  in  interpolating  the  yields,  allows  us  to  ascribe  a  thick¬ 
ness  to  the  response  surface.  To  illustrate  this,  the  zone  of  suitable 
performance  has  been  removed  from  the  cube  in  Photo  4  (Appendix). 
The  zone  of  suitable  performance  now  appears  as  three  boards  nailed 
together.  The  hollow  space  is  known  as  the  zone  of  acceptable  vari¬ 
ables  . 

Performance  confidence  requirements,  reliability  requirements, 
and  the  experimental  data  determine  the  thickness  of  the  surfaces. 
Only  one  way  of  applying  this  method  is  illustrated.  The  response 
surface  of  finite  thickness  would  be  used  to  construct  the  zones  in  dif¬ 
ferent  ways  for  different  performance  requirements.  Suppose  the 
velocity  were  required  to  be  108  ±  5  fps  instead  of  simply  greater  than 
108  fps,  still  keeping  the  pressure  and  acceleration  requirements  as 
before.  Then  the  zone  of  suitable  response  would  be  represented  by 
the  green  board  marked  V  in  Photo  4,  The  zone  of  acceptable  vari¬ 
ables  would  be  represented  by  a  surface  running  along  the  board, 
bisecting  the  thickness.  There  is  an  extremely  wide  variety  of  re¬ 
quirements  that  can  be  treated  with  this  technique.  No  unusual  or 
exotic  statistical  ma,thematics  is  required. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Our  general  conclusion  is  that  the  use  of  factorial  type  experi¬ 
mental  design  programs  represents  a  definite  advantage  to  the 
ballistic  designer.  These  advantages  are  measured  in  terms  of  a 
larger  number  of  variables  investigated  for  fewer  rounds  (time  and 
money  economy).  In  addition,  interaction  effects  among  the  vari¬ 
ables  are  determined.  Adding  the  Box  technique  and  pictorial 
representation  to  the  use  of  factorial  experiments  in  ballistic  re¬ 
search  gives  the  experimenter  a  more  economical  and  vivid  picture 
of  how  the  variables  operate.  To  this  picture  may  be  added  the 
variances  of  each  response.  Thus,  a  zone  of  suitable  performance 
may  be  determined  in  which  the  greatest  reliability  of  operation  is 
obtained. 

In  the  design  of  ballistic  devices  it  is  recommended  that  factorial 
experiments  be  conducted  and  combined  with  a  "box"  representation  of 
the  results. 
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